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Abstract 

The XXZ quantum spin chain has a triple point in its ground state h-l/A phase 
diagram. This first order critical point is located at the joint end point of the two 
second order phase transition lines marking the transition from the gapless phase to 
the fully polarized phase and to the Neel ordered phase, respectively. We explore 
the magnetization and the short-range correlation functions in its vicinity using the 
exact solution of the model. In the critical regime above the triple point we observe a 
strong variation of all physical quantities on a low energy scale of order 1/ A induced 
by the transversal quantum fluctuations. We interpret this phenomenon starting from 
a strong-coupling perturbation theory about the highly degenerate ground state of the 
Ising chain at the triple point. From the perturbation theory we identify the relevant 
scaling of the magnetic field and of the temperature. Applying the scaling to the exact 
solutions we obtain explicit formulae for the magnetization and short-range correlation 
functions at low temperatures. 

1 Introduction 

Critical points in the ground state phase diagram are nowadays often called quantum 
critical. In this terminology the triple point of the XXZ chain is a quantum critical point 
of first order. Quantum critical points have been proposed as a paradigm in condensed 
matter physics for explaining severe deviations in the behaviour of strongly correlated 
electrons from the Fermi liquid picture. Much of the intuitive understanding of quantum 
critical points comes from simple exactly solved models In this work we would like 
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to extend the list of those quantum critical points that can be studied exactly by another 
example. 

The ground state phase diagram of the XXZ chain in the h-l/A plane has three distin- 
guished critical points. Two of them are the isotropic ferromagnetic and antiferromagnetic 
points at h = and A = — lorA = l, respectively. Although the model is exactly solvable 
by Bethe ansatz, it is hard to calculate its thermodynamic properties in the vicinity of the 
ferromagnetic point. The antiferromagnetic point is easier to access, but still requires a 
certain effort in the numerical solution of the integral equations [13] . 

A third critical point, to which less attention has been paid, is the aforementioned 
triple point. It will be studied in this work. The triple point is located at A = oo and 
h = h c with A denoting the anisotropy parameter of the XXZ model and h c the critical 
magnetic field separating between the fully polarized phase and the Neel phase of the 
Ising chain. This critical point is simpler than the other two, because degenerate strong- 
coupling perturbation theory can be applied in its vicinity. Applying the same perturbative 
expansion to the integral equations that determine the thermodynamics of the model [13] 
we obtain free fermion type equations for low temperatures and algebraic equations in 
the zero temperature limit. Those equations provide explicit results e.g. for the ground 
state magnetization and some static correlation functions up to first order in 1/A. On the 
other hand, treating the integral equations numerically, we obtain the exact temperature 
dependent picture, and we can assess the range of validity of the perturbation theory. 



2 XXZ chain in the vicinity of the strong coupling critical 
point 

We consider the XXZ Hamiltonian, which on the 'gapped antiferromagnetic side' of the 
phase diagram is most naturally written as 



L r j h 



(1) 



with J, h > and A > 1. Here the aj, j = 1,... ,L, a = x,y,z, act locally as Pauli 
matrices, J is the exchange coupling, h denotes the strength of the magnetic field, and 
A = ch(r/) is the anisotropy parameter. 

The ground state phase diagram in the h-l/A plane can be obtained from the Bethe 
ansatz solution of the model (see e.g. [12]). We show a detail of it in figure [TJ The lower 
critical field h e (left red line) is given by the excitation gap [6j[T3] whereas the saturation 
field (right red line) is given by the simple explicit formula h s = 4J(1 + 1/A). These fields 
determine the second order phase transition lines between the critical and massive phases 
drawn in figure [TJ The lines are asymptotically straight. For the lower critical field one 
can see deviations from the asymptotic behaviour h e « 4J(1 — 2/A), but they are small 
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Figure 1: (Color online) Ground state phase diagram of the XXZ chain in the h-l/A plane 
for fixed J and large A in the vicinity of the triple point. 

on the scale used in the figure. The point h = AJ, A = oo in which the two second order 
phase transition lines join is the triple point we are interested in. In the following we shall 
characterize it starting with the Hamiltonian. 

For A — > oo the XXZ chain (JTJ) simplifies to the Ising chain 

1 L 

= g X, dia s( _/t > ~ 4J ' ~ 4J ' h )i-^i ' ( 2 ) 
j'=i 

which is diagonal in a basis of tensor products of local S^-eigenstates. From the second 
equation we can easily understand how its ground states change with the magnetic field. 
If h < AJ the two Neel states consisting of alternating up- and down-spins are the ground 
states. If h > 4 J the ground state is the unique fully polarized state with all spins pointing 
upwards. In the former case the ground state magnetization per lattice site vanishes, and 
in the latter case it is one half. Thus, there is a first order phase transition point at 
h = h c = 4 J. 

What happens exactly at the phase transition point? As can be seen from ([2]) all states 
with no neighbouring down-spins are degenerate and have lowest energy —2JL. Clearly 
their number is macroscopic in the thermodynamic limit where it is measured as the average 
number s of states per site. We can obtain it from the known [3] entropy per site in the 
zero temperature limit. It is equal to the golden mean, s = (1 + y/5)/2. All states with 
exactly two neighbouring down spins are the lowest excited states above the ground state 
with excitation energy AJ. The magnetization per lattice site right at the critical point 
takes a value different from zero and one half, namely m = l/2\/5. 
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Let us rewrite the Hamiltonian ([I]) as 

L 



J 

3=1 



with a = A(h/h c — 1). Then, for every fixed a, the second term becomes small for 
sufficiently large A and can be viewed as a perturbation to the Hamiltonian !K/(/t c ). The 
splitting of the ground state under the influence of the perturbation can be understood 
by means of an effective Hamiltonian obtained in second order degenerate perturbation 
theory, in a similar way as e.g. the Heisenberg Hamiltonian is obtained from the Hubbard 
Hamiltonian at half band-filling in the strong coupling limit (see e.g. appendix 2. A of [7]). 
Note that curves of fixed a are straight lines originating from the critical point at h = h t 
and A = oo in figure [Q 

Using second order degenerate perturbation theory we obtain the effective Hamiltonian 

X = -2JL + -P ]T (a^a* + 0)^0) - 2aaf) P 
i=i 

J L 



c 



j=l 



Po, (4) 



describing the splitting of the ground state. Here Po is the projector onto the subspace 
with no neighbouring down-spins. Note that the magnetic field does not enter into the 
second order correction. 

The first order term, let us call it "K\, is an integrable Hamiltonian [2] in full analogy 
with the situation with the one-dimensional Hubbard model in the strong coupling limit. 
The operator p\ = e --^i/ATy ^ T e -JJCi/AT re pi aces th e statistical operator to lowest order 
in J I A. For T — > 0+ it becomes the projector onto the ground state of "K\. It follows that, 
to leading order, the ground state correlation functions depend only on a, i.e. they are 
constant on straight lines originating from the triple point. Their values on lines parallel to 
the /i-axis in figure [1] must vary continuously between the values taken at the boundaries 
between the massive and the massless phases marked by the second order phase transition 
lines. The value of the magnetization per lattice site, in particular, calculated in the limit 
A — > 00 on straight lines of fixed a, must vary continuously between zero and one half, 
which already gives us a precise picture of the nature of the singularity of the first order 
quantum critical point. 

Considering the finite temperature physics of the effective Hamiltonian ^ for fixed a 
we see that, to leading order, the anisotropy A appears only as a prefactor, i.e. as a new 
temperature scale. This can be utilized by introducing a rescaled temperature r = AT. 
Keeping a and r fixed and sending A to infinity removes all corrections to "K\ . In this limit 
the original XXZ Hamiltonian is replaced by IH4. The form of the effective Hamiltonian 
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Figure 2: Magnetization of the XXZ chain close to the triple point for T — )■ 0. 

implies that this should be justified as long as A is large enough, say A > 10, since then 
the second term in becomes small compared to "Ki. At the same time the temperature 
T must be small compared to the excitation gap 4 J of IK/(/i c ) for the perturbation theory 
to be applicable. 

As we have seen, our strong-coupling analysis provides a complete qualitative picture of 
the ground state correlation functions in the region of the phase diagram above the triple 
point. Clearly, their values in leading order can be calculated using the exact solution 
of "Ki. Here we will proceed differently. We will take our previous results for the XXZ 
chain [13], already conveniently formulated in terms of integral equations, and consider 
them in the zero temperature limit and in the scaling limit. This way we obtain new and 
explicit results for the magnetization and some short-range correlation functions. They are 
shown in the following sections. 

The magnetization and the short-range correlation function in the full phase diagram 
of the antiferromagnetic XXZ chain (A > —1, h and T arbitrary) in the thermodynamic 
limit can be obtained by means of the formulae derived in [HE]. This only requires to 
solve certain well-behaved linear and non-linear integral equations numerically and can 
be done to arbitrary numerical precision. Examples were worked out in [U|13]. In |13j . 
in particular, we noticed the extreme variation of the short-distance correlation functions 
close to the triple point. Meanwhile we have worked this out in more detail. 

In [13] we discussed the temperature behaviour of two-point correlation functions, e.g. 
of the connected correlation function (cfcrg) — (erf) (03). It turns out, however, that the 
signature of the quantum critical point is most clearly observed in the magnetization per 
site (which is the only non- vanishing one-point function) , as it is monotonic in the magnetic 
field. We show it in figures [3] and [2l In addition, the two-point functions for two and three 
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Figure 3: (Color online) The one-point function (a z ) for large values of the anisotropy, 
fixed J and different values of h/J. The labels in the panels are the tuples A, h/J where 
'Ising' denotes the analytic Ising curves. 

sites are shown in figure |4l Note that the non- vanishing of the transversal correlation func- 
tions in the zero temperature limit is due to the residual quantum mechanical interactions 
for finite A. 

The left panel in figure [3] shows how small changes in the magnetic field or in the 
coupling 1/A induce large changes in the low temperature behaviour of the one-point 
function (cr|) due to the proximity of the triple point. The critical cone above the triple 
point corresponds to a G [—2, 1] (see below). The point A = 1000, h = 4.1 is outside the 
critical cone on the massive fully polarized side, while A = 1000, h = 3.9 corresponds to 
a point on the massive antiferromagnetic side of the phase diagram in figure [TJ All other 
curves in the left panel of figure [3] belong to various values of a £ [—2, 1]. The curves for 
h/J = 4, i.e. for h = h c , in particular, correspond to a = 0. It can be seen in the figure 
how, for increasing A, these curves become more and more similar to the curve for the 
Ising model at the critical point, but for small enough temperature always deviate from 
the value (crj) = 1/V~5 of the one-point function of the Ising chain. 

In the right panel of figure [3] we compare curves with different values of the magnetic 
field and a = 0.5 and curves for the Ising model with the same values of the magnetic field. 
We can identify four different temperature regimes. For very high temperatures T 3> h the 
curves are independent of A and h. Then, for intermediate temperatures T ~ h, they are 
independent of the (large) anisotropy and only depend on the magnetic field. This is where 
the curves for finite A and the curves of the Ising model match. We expect that this regime 
extends down to T ~ 10J/A, where the thermal fluctuations still dominate the transversal 
quantum fluctuations. Next comes a regime T ~ <//A, where thermal fluctuations and 
transversal quantum fluctuations are of the same order of magnitude. Here the correlation 
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functions depend on the anisotropy and the magnetic field. The different curves in this 
regime match, if one introduces the rescaled temperature r = TA, see section HJ Finally, 
for very low temperatures T <C J/ A, the product a = A(h/h c — 1) determines the value 
of the correlation functions. This case will be worked out explicitly in the next section. 
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Figure 4: (Color online) The temperature dependence of the two-point functions for two and 
three neighbouring sites are shown for different values of the magnetic field for A = 1000. 
Additionally, the zero temperature asymptotic values obtained in (|15p are shown as black 
lines. 
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3 Zero temperature limit 

In the following we shall work out explicit formulae for the zero temperature asymptotics 
of the magnetization and of a few neighbour two-point functions on lines of constant a = 
A(h/h c — 1) in the vicinity of the triple point. We shall resort to our previous work |13j 
in which we studied the XXZ chain for A > 1. The simple idea is to perform the zero 
temperature limit for large A and fixed a in the special functions of [13] that characterize 
the correlation functions. For this purpose the form of the non-linear integral equations as 
used in [5] is slightly more convenient than that of [13], for, in the zero temperature limit, 
where the integral equations become linear, the range of integration is over the Fermi sea 
and not over its complement. For large A, finally, the linear integral equations turn into 
algebraic equations. 

Referring to j8j,[T3] , where certain auxiliary functions b and a were explained, we in- 
troduce the familiar dressed energy e obtained in our setting from b or a in the zero 
temperature limit, 



e(x) = lim Tlnb(x) = lim — T mafia; — 

T->0 T->0 



(5) 



It satisfies the integral equation (see e.g. [8]) 



A 



e(x) = h-AJsh(r])K v/2 {x)- / — K v (x - y)e(y) 




(6) 



-A 



with 



K n (x) 



sh(2r?) 



(7) 



2 sin(x + \rf) sin(x — \rj) 



and e(±A) = 0. 

Now we use h = h c (l + t ) and expand the equation up to the order 1/A. Then 



A 




(8) 



which can be solved explicitly, 




(9) 



with c = ^A. The Fermi points ±A or ±c, respectively, are determined by the rescaled 
magnetic field a, 

a = (1 + c) cos(7rc) sin(-7rc) . (10) 
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Here c is a monotonic function of a for a G [—2,1]. The boundaries —2 and 1 of the 
interval correspond to the phase transition lines (red lines in figure [T]) determined to the 
order 1/A. For magnetic fields smaller than the lower critical field h e or larger than the 
saturation field h s the ground state and, hence, the Fermi points are independent of h. 
The ground state energy e is given by 

A 

h , f dy 



+ / ^ Kr]/2 {y)e(y) (11) 



7T 



-A 



h c f a\ h c ^sm(irc) - ac 

2 \ + A) A 1 + c • 1 j 



Using it one obtains for the magnetization in leading order 



(a*) 1-c 

m = ± — '- = — . 13) 

2 2(1 + c) K ' 

Then the leading order magnetic susceptibility is 

^ /i c (1 + c) 3 sin(-7rc) ^ ^ 

We see that it diverges at the critical point. 

In a similar way one can take the zero temperature limit for the functions that determine 
the correlation functions, e.g. for uj and u/ defined in [13] and for the auxiliary functions 
stemming from linear integral equations (see the appendix, where the zero temperature 
limit is obtained from the low temperature approximation discussed in the next section). 
This leads to the following expressions for the zero temperature correlation functions 

1 — V 
2sin(7rc) 

i^) h = -^ T ^, (15b) 
. z l-3c + 4c 2 4sin 2 (vrc) , 

=^T-c ^TTcT ' (15c) 

2sin 2 (7rc) 2csin(27rc) sin(27rc) 
( " lCJ3)ft = ^OTcT " vr(l + c) + ~^r~ ■ 

In figure|3]we show that the asymptotic values of the zero-temperature correlation functions 
match indeed the low temperature behaviour of the temperature dependent correlation 
functions calculated numerically from the integral equations as described in |13j . 

In figure [5] we show the dependence of the correlation functions as calculated in (|15p 
on the rescaled magnetic field a. In addition the normalized Fermi point c is shown as 
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Figure 5: (Color online) The asymptotic values of the zero temperature correlation func- 
tions according to equation (|15h are shown as a functions of the rescaled magnetic field a. 
In addition we show the normalized Fermi point c as a function of a. 



a function of a. Both, the Fermi point c and the one-point function (cr z ), are monotonic 
functions of the magnetic field. The connected two-point function (afa^) — (erf) (<r|) is 
monotonic as well. The other correlation functions in the figure, however, show non- 
monotonic behaviour. Both two-point functions have their extremum above the critical 
field of the Ising chain, whereas the minimum of {crfa^) is exactly at the critical field of 
the Ising chain. 

4 Low temperature behaviour 

After the discussion of the zero temperature limit in the previous section which gave the 
correct scaling behaviour of the magnetic field, we would like to discuss small but finite 
temperatures. For this purpose we introduce a new rescaled temperature r = AT and 
expand the integral equations to lowest non- vanishing order in 1 /A with r and a kept fix. 
This leads to 

tt/2 

e{x) = h c a - h c cos(2x) + J ^ r In (l + e - £(y)/T ) (16) 

-tt/2 

where e{x) is the finite temperature generalization of ([5]) in the scaling limit. Therefore we 
use the same notation for this function. We see that only one interval [— 7r/2,7r/2] of the 
original contour in the non-linear integral equation for the full XXZ chain [5] survives the 
limit. The contribution of the remaining part of the contour is exponentially suppressed 
for A — > oo. 
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The solution of the equation (|16p is nothing but the dispersion relation of free fermions 
on the lattice 

e(x) = (J, — h c cos(2x) , (17) 

however, with a chemical potential \i which depends non-trivially on magnetic field and 
temperature and which is implicitly determined by the equation 



tt/2 



H = h c a+ J ^rln(^l + e- ( ^ cCOs(2y))/T ) 

-7T/2 

The free energy in this approximation is then given by 

vr/2 

/ = _^(l + £)+!. J ^ T l n (l +e -(^cos( 23/ ))/^ (1Q) 

as for free fermions. Defining 



-tt/2 



tt/2 

f dy cos{2ny) 

J TT l + e (M-^cos(2j / ))/r ' ^> 



-tt/2 



the magnetization follows as 



The calculation of the two-point functions within the same approximation is sketched 
in the appendix. We obtain 

1 + Jo 

(*M) T , h = - yf^ , (22b) 

l-3J + 4/ 2 47? 
<"i"3> Tlfc = 1 + Jo "j^, (22c) 

1 + Jo i + J o 

An important remark is that the structure of the two-point functions (|22|) . unlike the 
structure of the free energy, is not the same as for free fermions. They rather depend on 
the density-density two-point functions I n for free fermions in an interesting way. 
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Figure 6: (Color online) The longitudinal correlation functions for short distances as a 
function of the temperature are shown for different values of the anisotropy and a = 0.5. 



In figure [6] and [7] the short-range longitudinal correlation functions are shown for fixed a 
and different values of the anisotropy. The exact solutions are compared with the solutions 
obtained within the scaling approximation. The agreement is surprisingly good. The ap- 
proximate formulae hold up to T ~ J. This is compatible with the fact that the degenerate 
perturbation theory should be valid as long as the temperature is small compared to the 
excitation gap 4J. Alternatively we can understand this fact from the perturbative treat- 
ment of the integral equations. As stated above and explained in some more detail in the 
appendix, the contribution of the second contour to the integral equation is exponentially 
suppressed with A for small enough r. On the scale of the true temperature T this means 
that it behaves like exp(— J/T), i.e., in particular, in a way independent of A. 

In order to understand the strange high temperature behaviour of the correlation func- 
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Figure 7: (Color online) The longitudinal correlation functions for short distances as a 
function of the temperature are shown for different values of the anisotropy and a = —1.5. 



tions in the scaling limit, we have to recall that they are the exact correlation functions of 
the Hamiltonian "K\ , which is defined on the restricted Hilbert space with all configurations 
containing neighbouring down-spins excluded. In the high temperature limit, where statis- 
tical physics means to simply average over all possible configurations, the excess of up-spin 
electrons becomes visible in a non-vanishing magnetization. We can see in the figure that 
its high temperature value becomes independent of the anisotropy, as expected. Similarly, 
projecting out configurations with neighbouring down-spins means that it becomes more 
likely to have antiparallel neighbours, which causes the longitudinal neighbour-correlation 
function to be negative. 

The scaling approximation applies to a wide range of the anisotropy. For A = 100 and 
A = 1000 the deviation from the exact solution is hardly recognizable, and the approxima- 
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tion is still in good qualitative agreement with the exact solution for A = 10. Comparing 
the cases a = 0.5 in figure and a = —1.5 in figure El one sees that for A = 10 the offset 
in the zero temperature limit is larger in the latter case. This is in agreement with the 
observation that the critical lines in the zero-temperature phase diagram in figure Q] relate 
differently to a. The phase transition line to the fully polarized state belongs exactly to 
a = 1, while the line marking the transition to the Neel phase belongs to a = — 2 only 
asymptotically and deviations become visible for A ~ 10. For this reason we generally 
expect larger deviations for smaller a E [—2, 1]. 

In general both figures show that in the regime T ~ J/ A, in which the correlation func- 
tions depend on the anisotropy and the magnetic field, the dependence on the anisotropy 
can be removed by introducing the rescaled temperature r. Since the scaling approxima- 
tion is valid for temperatures up to T ~ J and since the XXZ chain is well approximated 
by the Ising chain for temperatures higher then T ~ 10 J/A, the critical regime of the XXZ 
chain above the triple point is fully accessible by simple approximations (either Ising or 
'Ki) as long as A > 10. Outside the critical cone the Ising approximation works well for 
A > 10 and all temperatures anyway. 

5 Conclusion 

Based on the exact solution for the thermodynamics and short-distance correlation func- 
tions |13j we have examined the vicinity of the triple point of the XXZ chain in the critical 
regime. The low temperature physics in this region of the phase diagram and the nature of 
the singularity at the critical point can be understood qualitatively be means of degenerate 
perturbation theory around the ground state of the critical Hamiltonian "Ki(h c ). From the 
form of the second order effective Hamiltonian (j3J) it follows that a small scale T ~ J / A ex- 
ists on which the correlation functions show strong variations and that below this scale the 
correlation functions must be constant on straight lines of fixed a = A(h/h c — 1). Using the 
exact solution, formulated in terms of certain linear and non-linear integral equations [13], 
we have obtained a full quantitative understanding of the vicinity of the critical point at 
any temperature. We have seen that, as long as A > 10 the full XXZ Hamiltonian in the 
critical phase above the triple point is well approximated by the effective Hamiltonian "Hi 
as long as T < J and by Hj(h c ) for T > 10 J/ A. In particular, for the interesting interme- 
diate and low temperature regime T < J it is enough to consider "K\. This is interesting, 
since "K\ belongs to the relatively simple class of Hamiltonians similar to the impenetrable 
Bose gas [2]. The Bethe ansatz solution of models in this class is simple enough to admit 
for the derivation of finite temperature dynamical correlation functions in a closed form 
called determinant representation [ID]. A determinant representation for "K\ was derived 
in [TJ. It could be the future starting point for an analysis of the asymptotic behaviour of 
the dynamical correlation functions in temperature regions in which conformal field theory 
is no longer a valid approximation. Such an analysis would be complementary to the per- 
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turbative study of the finite temperature dynamical correlation functions in the massive 
Neel ordered regime recently presented in [5]. 
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Appendix: Low temperature approximation and zero tem- 
perature limit of the linear integral equations 

We are referring to the notation and definitions of |13J . In order to calculate the low 
temperature limit of the two-point functions we have to begin with the functions G and 
G' . For these functions a shift as in ([5]) has to be made. We define 

9u W : = -G(ix - rj/2, w), g „ (x) := (23) 

which is consistent with the notation used in |13j except for the denominator r/ in the 
second equation introduced for convenience here. 

In the low temperature limit for large A the linear integral equations for G and G' 
imply 

tt/2 

g-(x) = -2K v/2 (x -„)-/ ^K v {x - y ) —^M—- , ( 2 4a) 

-7T/2 

tt/2 

f d V t t \ 9u(y) 



where 



1 _|_ e (jM-h c ooa(2y))/T 
tt/2 

tt/2 

—KJx - y) 7^^- , (24b) 

7T ,?v yJ \ _(_ e (At-h c cos(2y))/r ' > 

tt/2 



L v (x) = , iSi " (2x) , r . (25) 

2 sm(x + 177) sm(x — lrj) 



The second part of the integration contour vanishes in this limit as it gives only exponen- 
tially small contributions for A — > 00, similar to (116p . 
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We solve these equations up to the order exp(— 2rf) and obtain 



9u ( x ) 



1 + it 



+ e" 



4/i 
l + /t 



cos(2zv) — 4cos(2(x — u)) 



+ e" 



-2'/ 



4J 2 4J? 4Ii 

— — - cos(4z/) - 2 - 4 cos(4(x - 1/)) + — — cos(2x) 

1 + J (1 + Jo) 1 + / o 



(26a) 



1 + Jr 



+ e" 



2c 



-i2(a:-!/) 



2J 



1 + Jo 



1 j2u 



+ e" 



-2,, 



2c 



-i4(x— v) 



i4Ji , , 2Ji , , 2J 2 u „ 2J? 
+^^- sin(2x) - - — — cos(2x) - ^^e l4ly + 



(26b) 



1 + J ~"~ V ~ ' 1 + J ~~~ V ~ ' 1 + J ~ ' (1 + Jo) 2 . 
Finally we need the low temperature limit of u> and ui' defined in |13j . One easily finds that 

tt/2 



r dx 

V2) = K n {v\ -D 2 ) + 2 I —K v/2 (x - u 2 )j 
^2) 



-tt/2 
tt/2 



-L v (i>i - z> 2 ) 



d.r 



_|_ e (n-h c cos(2x))/t ' 



(27a) 



-tt/2 



1 Ctg " P2 " ^ 1 + e (M-/W2*))/ T 
tt/2 



+ 2 J f^ 2 (x-^) 1+ /^ w)/T , (27b) 

-7T/2 

where we used the notation = — ii/j. Inserting (|26p here and using equations (24) and 
(25) of [13] we obtain the equations ([22]) . 

For the zero temperature limit integrals over the Fermi weight of the form 



tt/2 



dx 



7T I -L e (/-t—hc cos(2x))/t 



(28) 



-tt/2 



have to be changed into integrals over the Fermi sea 

A 



dx 



7T 



fix). 



(29) 



-A 



The functions J n , for instance, that determine the correlation functions in the low temper- 
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ature scaling approximation have the following zero temperature limit 

Io — > c , (30a) 

h -+ ^ , (30b) 
sin(27rc) 

h -> — ^ • (30c) 

With this limit introduced into (|22p one obtains the zero temperature correlation functions 
(El). 
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